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ABSTRACT 

Quantum dissipation in thermal environment is investigated, using the path inte- 
gral approach. The reduced density matrix of the harmonic oscillator system coupled 
to thermal bath of oscillators is derived for arbitrary spectrum of bath oscillators. 
Time evolution and the end point of two-body decay of unstable particles is then 
elucidated: After early transient times unstable particles undergo the exponential 
decay, followed by the power law decay and finally ending in a mixed state of resid- 
ual particles containing contributions from both on and off the mass shell, whose 
abundance does not suffer from the Boltzmann suppression. 
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It is a long standing problem in physics to resolve how a quantum system evolves 
with time under influence of an environment described by some mixed state. The 
effect of the environment on the system dynamics is governed by the fluctuation 
(or the noise) and the dissipation, related by the fluctuation-dissipation theorem. 
The simplest, yet the most fundamental model of quantum dissipation is harmonic 
oscillator coupled to infinitely many oscillators that make up a bath in a mixed state. 
Two powerful methods to analyze this problem are the quantum Langevin equation 
|TJ], and the path integral approach M, ||, M. 

The important and the difficult part of analysis of this problem is how to treat 
the non-local part of the correlation of the system variable obtained after integrating 
out environment variables. The non-locality is required from quantum mechanical 
principles. Many past studies have relied on the local form of the dissipation kernel in 
the path integral approach. Some exceptions are consideration in limited applications 
such as || , H . It has been pointed out that this often-used approximation of the local 
friction (in the path integral sense) has a fundamental flaw; violation of the positivity 
of the density matrix J7|. This has made it difficult to extract in a reliable way a 
fully quantum mechanical time evolution of the small system immersed in thermal 
environment at early times, especially when the temperature is low. Needless to 
say, the approximation of the local friction is excellent in the most important phase 
of time evolution, thus is useful in many practical applications. But as we shall 
see, there are physical effects that can be understood only with the full non-locality 
incorporated. 

In the present work we extend previous works on the non-local dissipation Q 
and give a quantum mechanical formulation of the system evolution in thermal en- 
vironment. The present model itself is not new, and indeed, is a standard model 
of system-environment interaction discussed repeatedly in the quantum dissipation 
problem. A new method is developed by fully exploiting analyticity familiar in scat- 
tering theory, to treat the general class of environment models characterized by the 
spectral weight consistent with general principles, but otherwise taken arbitrary. We 
then apply this formalism to the two-body decay of unstable particle in field theory, 
by identifying as the environment oscillator a composite operator of two decaying 
particles. We may thus discuss details of the decay law and what is left behind after 
decay. It is the extended analyticity of a kernel (denoted later by F(z)) that allows 
us to treat the whole time region of the decay process in a unified way. 
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Another difference of our study from most of the past ones is our assumption 
of the presence of a gap in the spectrum of environment oscillators. Although this 
appears a minor technical point, it actually leads to differences at low temperatures 
and also for behaviors of physical quantities at very late times. The gapless case 
often appears in condensed matter physics, for instance phonons in medium. But 
there are cases in which one cannot neglect the gap such as in the unstable particle 
decay if the mass of the daughter particle is finite. 

Let the system variable in question be denoted by q and the environment variable 
by Q a - The Lagrangian of our problem consists of three parts: 

L = L q [q] + L int [q,Q} + L Q [Q}. (1) 

We take for the system-environment interaction the bilinear term with infinitely 
many harmonic oscillators for the environment: 

Lint = - 9 S Ca Q a ' ( 2 ) 
a 

L q = ^(q 2 -u; 2 q 2 ), (3) 

L Q = W(Ql-^QD- (4) 

Z a 

The renormalization of the frequency will be discussed later in appropriate places. 
Introduction of a counter term in relation to the renormalization as often discussed 
in the literature is actually the problem of how to relate the parameters of the 
theory such as uj to observable quantities. We shall discuss the relevant observable 
parameter Hi later. 

The response of the environment to the system is characterized by what we 
call the response weight, denoted by r{uj): For the dynamical system of harmonic 
oscillator of frequency ui coupled to environment oscillators of discrete frequency 
u a , it is 

r(ui) = £ |j- S(ui - ui a ) = jT dui *M D(ui) , (5) 

where in the continuum limit D{uS) is the density of states per unit frequency. ui c 
is the threshold for r(uo) ^ 0. Although we often write expressions in the discrete 
distribution of frequencies, we have to consider the continuum limit. It is often useful 
to extend r(u) to the region ui < by r(— ui) = —r{uj) . 
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It is very useful to introduce the influence functional J-"[g(r) ,q'(r)], following, 
but slightly modifying, the original definition of Feynman and Vernon JJ]. The 
basic idea is that one is interested in the behavior of the q— system alone and traces 
out the environment variable altogether. This way one can compute the reduced 
density matrix for the system that contains all informations without knowing too 
much details of the environemnt. We define the influence functional by convoluting 
with the initial state of the environment unlike the original one. Once the influence 
functional is known, one may compute the transition probability and any physical 
quantities of the q— system by convoluting dynamics of the system under study. For 
instance, the transition probability is given, with introduction of the density matrix 
P {R) , by 

J Vq(r) J Vq\r) J dq { J dq[ J dq f J dq' f 

nq f ) nq'i) FU{t) , q'(r) ] e iS ^ iS ^ ^(q' f ) 

= Jdq f J dq' f r(qf)^(q' f )p iR) (q f ,q' f ), (6) 

where t/Vs are wave functions of the initial and the final q— states, and S[q] is the 
action of the q— system. 

The form of the influence functional is dictated by general principles such as prob- 
ability conservation and causality. Feynman and Vernon found a closed quadratic 
form consistent with these, 



nq(r),q'i 



T 



exp 



- jf* dr £ ds ( ^r)a R (r - s)£(s) + % £(r)a,(T - a)X(s) 
with £(r) = q(r) - q\r) , X{r) = q(r) + q'(r) . 



(7) 



Thus two real functions c^tV) are all we need to characterize the system-environment 
interaction. These are defined here in the range of r > . The fact that «j depends 
on the difference of time variables, r — s, is due to the assumed stationarity of the 
environment. 

The correlation kernels appear in the influence functional as a form of the nonlocal 
interaction and they are the dissipation aj and the noise a^, related to the response 
weight by 

ot/(r) = - - / dcur(to)e- iujT , (9) 



2 

a R (r) = \ r dwcoth(^)r( U )e-- . (10) 
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In the path integral approach integration over the sum variable X(t) is trivial 
for the harmonic oscillator system, since both the local part and the nonlocal action 
above are linear in this variable: 

111 (^)X(r)-c 2 e(r)X(r))- 

f dr f T ds ( Z{r)a R {T ~ s)£{s) + i £(r)a/(r - s)X(s) ) . 
Jo Jo 

Thus result of the path integration of the system variable X(r) gives the classical 
integro-differential equation that £(r) must obey 

+ Wo 2 £(r) + 2 £ dst(s)ai(s-T) = 0. (11) 

The end result of the £ path integral then contains an integral of the form, 

- f dr [ T dst(T)a R (T-s)S(s), (12) 
Jo Jo 

using the classical solution £ (r) with specified boundary conditions, £(0) = & , £ (£) = 

£/■ 

In the local approximation often used the dissipation kernel of the form is taken 
to have the form of 

a/(r) = 5u 2 6(t)+t}5'(t), (13) 

with 5u 2 the frequency shift and the r] term representing the local friction. This 
choice enables one to solve the £ equation (|TT| ) by elementary means. On the other 
hand, the noise kernel is given by the response weight of the form, 

r(u,) = ^u,/(£), (14) 

with f(x) some cutoff function and Q a high frequency cutoff. At high temperatures 
this approximation reduces to the well known classical form of the fluctuation, 

Mr) = — S(t). (15) 

7T 

At low temperatures, howerver, the use of the cutoff function only in the noise kernel 
while retaining the local form of the dissipation makes validity of this approximation 
dubious. Our present approach does not make this local approximation, and instead 
uses exact solutions of the classical equation (|ITD . This equation has been used in 
other approaches ||, 0, |Q, too. 
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Two observations are crucial for our subsequent development: First, the integro- 
differential equation (|TT|) is converted to the one having the Volterra type kernel 
with the variable change r — *> t — r . It can then be solved by the Laplace transform, 
if the Laplace transform of aj(r) denoted by 6c(p) has a simple form. Second, in 
thermal environment a(r) = or(t) +zaj(r) is identified || to the real-time thermal 
Green's function, hence its Fourier transform has the well known analyticity property. 
The origin of the analyticity is ultimately traced to the causality: any physical 
disturbance cannot propagate faster than the light. 

Let us explain this in more detail. First, the Laplace transformed kernel of 
evolution operator is 

p 2 + co>g + 2a(p) 



with real p > 0. 



a(p)= 1 -r d u^± (17) 



is equal to the boundary value of the analytic function TrG(ip). This analytic function 
G(z) is related M to the imaginary-time thermal Green's function || and is regular in 
the complex z plane except on the real axis where it has the branch point singularity 
with the discontinuity given by r(uS). Our notational convention here is that uj and 
p are real and the complex variable is denoted by z. 

We next introduce the analytically extended function by 

FW = M - feH __L__ (18) 

TK>)=rT—- < i9) 

J-oo 111 Z — UJ 

It can be proved from the positivity of the response weight r(u) supplemented by 
the spectrum condition later explained that F~ 1 {z) has no complex zero. ||1 1|| Hence 
F(z) is regular except on the real axis. Its discontinuity across the real axis 

H(uj) = (F(uj + it) - F(uj - it) ) /(2m) (20) 

is fundamental to our subsequent analysis. In terms of two real functions, n(u;) and 
T(uj) defined by 

2nG(uj + ie) = U(u) -iuT(u) , (21) 

rr (l x _ I ^V) . 99 A 

1 ' ~ vr ( ^ - - n(w) f + to* T^cu) ■ [ ' 
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Furthermore, r(uj) = uoT(uo)/tt , and the real part IL(u) may be written by the well 
known dispersion integral: 

n(w) = V r duo' ^-L , (23) 



with V denoting the principal value of integration. 

As in scattering theory ||10|| , the analytic function F(z) may be extended to the 
second Riemann sheet through the discontinuity formula, 

F(uo + ie) - F(uo - ie) = i2<Kr(uj)F(u + ie)F(uo - ie) . (24) 

In the second lower sheet one finds a finite number of poles at z obeying 

z 2 - uo 2 - 2tt G(z) + 2nir (z) = . (25) 

Except at these poles F(z) is regular in the unphysical sheet. 

Physical interpretation of these singularities is possible in terms of the spectrum 
of the entire harmonic system. Due to the mixing with the environment the system 
spectrum is modified and is given by an eigenvalue equation for A = uo 2 . When 
the response weight r(u) — > a constant at uo — > oo, as happens often and in field 
theory models later discussed in particular, one has to subtract a term from integrals 
containing the response weight. This corresponds to renormalization of the bare 
frequency ujq with the frequency shift, 

POO 

5uo 2 = -2 / duor{uo)/uo. (26) 

J u c 

The renormalized frequency is given by uo 2 R = uo 2 + 5uo 2 . The eigenvalue equation 
written using uo 2 R is then 

2 , f°° , 2r(uo) 



^-^t^r - < 27 > 



U> c 



The left hand side coincides with — F The condition of stability of the entire 

system requires that the smallest eigenvalue A > 0, giving uj r > , or equivalently 



uOq>2 duor(uo)/uo, (28) 

if the integral on the right hand side is convergent. This is the condition needed to 
ensure the analyticity of F(z) as described above. 
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We do not consider the possibility of an isolated pole of F(z) on the real axis 
outside the cut, which corresponds to a stable state, implying that in the problem 
raised here a prepared initial configuration of the system does not completely decay. 
We thus focus on physical systems that have a continuous spectrum alone when 
interaction between the small system and the large environment is considered. 

The complex pole in the second Riemann sheet describes the behavior of the 
system oscillator interacting with environment oscillators: The real part of the pole 
position uj 2 gives a physical frequency squared including the frequency shift, and its 
imaginary part is related to the decay rate of the system oscillator. The pole mass 
here uj ^ uj r . Note that the pole position uj is, but ujr is not, an observable quantity. 
We assume as a physical requirement that a single pole exists in the nearby second 
sheet, but it should be evident that an arbitrary number of poles can readily be 
incorporated. 

Subsequent computation frequently uses the Laplace inverted evolution operator 
g(r) which can be written in terms of H(u): 



DC 



g(r) = 2 / dLuH(uj)sm(LUT). (29) 
Solution of the integro-differential equation (O) is then given by 
fH e s(t-r) ( 9(t-T)g(t) \ 

with the dot denoting derivative. Both g(t) and g(t) can be shown to satisfy the 
integro-differential equation of the form (x = g or g), 



-^- + uj 2 x + 2 f dr otAt - t) x(t) =0. (30) 
dt z Jo 

The reduced density matrix of the quantum system at any time is obtained 
from the action written in terms of the boundary values, S c i{£f ,Xf ; £j , JQ) , by 
convoluting with the initial density matrix. This action is computed as 

is* = -^q-^g-wtitf+iXfif-^Xiii, (si) 

U = 2 f dr [ T ds z(r) a R {r - s) z{s) , (32) 



o 



V = 2 dr dsy(r)a R (r-s)y(s), (33) 
Jo Jo 

W= f dr f ds(y(r)z(s)+y(s)z(T))a R (T-s) } (34) 



8 



*) - ^ ^ (35) 

<r) = g(t - r) - g(t - r)^|, (36) 
fiM c 9(t ~ r) , g(t - r)g(t) \ 

ar) = -*-jyr-* , [ 9{t - T) — ^^J- (37) 

The same effective action as ours has been derived by Grabert et al., as summarized 
in H . Our derivation here simplifies their calculation, with a new form of the reduced 
density matrix, Eq. (j38"D that follows, by fully exploiting the extended analyticity. Our 
method also makes possible a unified treatment of the exponential and the power 
law decay, as will be made clear later. 

We take as the initial state a product of thermal states, a system of temperature 
To = and an environment of temperature T = 1/(3. We may take To = T when 
we apply to the decay process of excited level initially in thermal equilibrium. In 
the limit of T — > it describes the ground state of the system harmonic oscillator. 
(We have recently computed the reduced density matrix, starting from another pure 
state, the first excited level of the system ocsillator. As expected, the late time 
behavior in this case is the same as in the present case.) 

After a series of straightforward Gaussian integration we find the reduced density 
matrix of the form, 

pW(X f , £/) = 2^exp[ -AX} - Bi) + iCXfa } , (38) 

A=—, B=~(I 3 -$), C=—, (39) 
8ii ' 2 v 3 h ' ' 2Ji ' K ' 

h = I[\h(u ,t)\ 2 ] + ±- coth(^)(<) 2 + c o y) , 

h = m\ h(u , t)k*(u ,t) ] + coth(^)^ W g) , 

I 3 =l[\k(u,t)\ 2 ] + ^oth(^)(g 2 W Q g 2 ) . (40) 

uq is a reference frequency taken as the initial system state, and equated here to the 
initial system frequency. If one so desires, either the renormalized ujr or the pole uj 
may be taken as another choice. But we imagine the situation a small system was 
added to a large environment at some time, its mutual interaction being absent prior 
to the initial time. In this circumstance it is appropriate to take uq as the reference 
frequency. As will be made clear later, dependence on the initial state dies away 
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quickly as time passes. In writing the reduced density matrix, we introduced the 
following notations; 



X[ /(«) ] = j°° du coth(^)r(u;)/(u;) , (41) 



and 



/i(w,t) = /' drg(r)e-^, (42) 
k(w,t)= [ drg(T)e- iu}T = g(t)e~ iult + tuoh(uj ,t) . (43) 

JO 

The density matrix p( R > from which any physical quantity can be computed at any 
time has explicitly been given by the discontinuity, H{uS) or r(u). 

This density matrix p( R > is positive definite and behaves acceptably at early times 
unlike the one in the local approximation, as will be explained elsewhere [O. Our 
interest here is limited to the late time behavior. 

It is important to know the behavior of the kernel g(r), which is found by de- 
forming the contour of u integration into the sum of the pole contribution (at z = Zq 
with Qzq < 0) in the second sheet and the contribution parallel to the imaginary 
axis passing through z = lu c , both in the first (I) and in the second (II) sheet JT3|: 



g(r) = \ Ke 



+ 3 



7T 



poo 

/ dy e~ yT ( F 1 {uj c + iy) - F n (u c + iy) 
Jo 



(44) 



with K" 1 = z Q — ttG'^Zq) + iTrr'(zo) . 

As seen from this formula, the pole contribution given by the first term describes 
the exponential decay, which usually lasts during the most dominant phase of the 
decay period, while the rest of contribution gives the power law decay at very late 
times; oc t~ a ~ l , and well-behaved early time behavior. The power — a — 1 is dictated 
by the threshold behavior of the response weight, r(u) oc (u — u c ) a . a = 1/2 for the 
S-wave two-body decay of unstable particle |0|. a > — 1 is required as a consistency 



of this approach; convergence of the u integration. 

It is often claimed that the power law behavior is never observable, since by 
the time this term dominates the exponential decay essentially eliminates the initial 
population. This is perhaps so when the environment is at zero temperature. In 
order to check observability of the power law decay in medium, we shall estimate the 
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transition time t* from the exponential to the power law period at low, but finite 
temperatures. Let us examine a typical case by taking the form of the response 
weight, r{uj) = c (lo — u c ) a , with < a < 1 in the range of u c < uj < Q (Q ^> u> c ) 
and with Co 3> Max ( uj c , T ) . The behavior of g(t) in the power law period is 

2c , , cos(uj c t + 77 a ) , , „ 

g (t)«-r(a + l) ; +1 2 (45) 

By equating this to the expression for the same quantity g(t) in the exponential 
period, one obtains 

u ™ \ ln ( 2 cr ( a+ 3 i) 7 ^ ) ' (46) 

with 7 = — Qzq the decay rate. In the weak coupling limit the factor inside the 
logarithm is large (oc c~ 2a ~ 3 ), and by the time the initial population has decreased 
like 

At non-zero temperatures the response weight r(u;) may depend on the environment 
temperature, as will be made more explicit in our application to the unstable par- 
ticle decay, for instance in eq. (|56|) . Temperature dependence of the parameters in 
these formulas is thus needed to check observability of the power law decay at finite 
temperatures. 

Emergence of the power law term in the quantum Brownian motion has been 
noted in some specific models of Ohmic type, [HI, (TB|, p!7| |, but we find this behavior 
as a general property in the presence of the non-local dissipation. For instance, the 



power law behavior of the two point correlation function was noted in ref. |L5|, |L6 



for specific models. We have instead demonstrated the power law behavior for the 
kernel function g(t) by separating the non-pole contribution that is essential to the 
power law behavior. Hence we have shown this property in more generality. We 
would like to stress that the presence of the branch cut singularity of the function 
F^ 1 (z) in the z plane is very important to derive the power law decay. This cut is 
further related to the non-analytic property of the real function r{uj) at the threshold, 
including the special case of the gapless u c = 0. For instance, if one takes the Ohmic 
form of the response weight given by pole terms alone, such as 

r«(*) = ^ , (48) 

(lu 2 - ft 2 + ^) 2 + fi 2 7 2 
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then one does not obtain the power law decay for g(t). For the gap less case of a 
fractional power a the formula fl45|) is still valid, but the kernel g(t) does not exhibit 
the power law behavior for an odd integer a, as seen from (|45| ) with uj c = 0. In this 
case one can define a regular function r(uj) including u = 0, thus the gapless case 
of odd integer a is exceptional in the sense that only this case does not produce the 
power law behavior. The critical condition for the absence of the power law behavior 
of g(t) is then regularity of r(u) at the threshold. 

The asymptotic late time behavior of the reduced density matrix is determined 
by h{uj , oo) = / °° dr g(r)e~ luJT , which is shown to be equal to the boundary value, 
F{uj — ie) , thus giving 

r(oo)\h(uj,oo)\ 2 = H(u) , (49) 
r(u)\k(uj,oo)\ 2 = uj 2 H(uj). (50) 

In these computations the 
asymptotic late times 

A - 

B - 
C - 

Note that the dependence on the initial state of the system via the /3q = 1/T factor 
disappears, hence no memory effect of the initial state remains in the asymptotic 
final state. 

In the high temperature limit 

^~^f, (54) 

This result is almost equal to the pole contribution, igoring a minor difference be- 
tween ujr and to. Contribution during the period of the power law decay is given by 
the continuous integral and is numerically subdominant, suppressed by the factor 
1/T 2 relative to the one from the exponential period. 

On the other hand, at low temperatures the contribution from the threshold 
region, uo m u c , cannot be ignored, giving the dominant contribution to the power 
law period. 



analytic structure and only that is important. Thus, at 



l^£duJCoM^)H(u)^ , (51) 

- / du coth(^) u 2 H{lj), (52) 
2 J uj c 2 

0. (53) 
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Let us apply these considerations to the decay of unstable particle; p> — > x + X ■ 
We assume that the decay product x is a P ar t of thermal components that make up 
the environment. The parent particle p may or may not be in thermal equilibrium 
with the rest of medium. Since we focus on the late time behavior, the initial state 
dependence disappears. Unlike the unstable particle decay in vacuum some amount 
of parent particles are expected to be left behind, even much later than the decay 
lifetime, because in thermal environment even heavier parent particles can be created 
by energetic particles of smaller mass. If so, what is the fraction of the parent particle 
left behind? 

The interaction of the system field p and the environment field x is assumed to 
be described by a relativistic field theory of Yukawa interaction Lagrangian density, 
| ^PX 2 ? with jj, a coupling constant of mass dimension. We presume that interaction 
among created x particles themselves is weak enough. In this case it is possible to 
identify as the environment variable the two-body bilinear operator of Xi 

E c ^ = f J d 3 x X 2 (x)e-^, (55) 

with k the momentum of p particle. Due to the assumed homogeneity each k mode 
can be treated independently. The continuous label u in the environment variable 
Q(oS) is identified here to the internal configuration of two body states with a given 
total momentum k. This field theory model was introduced in ref |Q. 

The relevant response weight for p — > x + X has been given in ref || . We shall 
recapitulate the main point of that calculation. The important point is that to lowest 
non-trivial order of the coupling /i, r(u) is given by the imaginary part of the one- 
loop diagram of the self-energy of p at finite temperatures, p being off the mass shell: 
uj 2 — k 2 7^ the p mass 2 . The result for the x loop diagram at finite temperatures 



is well known ||14|| , ||, and physically consists of two parts; the process p «-> x + X 



in the region of uo > \Jk 2 + 4m 2 and the other process p + x X (forbidden 
when all particles are on the mass shell, but allowed in thermal environment) in 
< uj < k, where m is the daughter x mass. The response weight riiS) does not 



vanish for |u;| > y k 2 + 4m 2 and |a>| < \k\ from the kinematics of the decay and the 
inverse decay of particles off the mass shell, with the constraint of the momentum 



conservation. Thus a gap of the spectrum exists in k < \u>\ < y k 2 + 4m 2 . The finite 
non-vanishing mass of the daughter particle (m ^ 0) is important for the existence 
of the gap and for associated physical consequences that follow. 
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For u > \Jk 2 + 4m 2 the response weight is 



a 2 ( Am 2 2 1 — e _/3a;+ 

(56) 

with 

, k 



W± = -±- v /l-4mV(^_p ) . (57) 

For < c<j < k only the second term in the bracket of Eq. (|56"D contributes. 

A useful, and adequate approximation we exploit for subsequent estimate is the 
weak coupling scheme with correct threshold and asymptotic behaviors incorporated: 

F{z) = I , . (58) 

In this approximation we replaced the real part TL(u) by the constant pole location 

U). 

A quantity of physical interest is the fraction of remaining particles given by the 
occupation number at asymptotic late times, 

1 / Pk . 2\ _ Bk . UJ k 1 



for each k mode = y k 2 + M 2 ). The temperature dependent part of this quantity 
is 2 

< = r du "^r — t + — ) H ( u ' fc ) > ( 6 °) 

Jo e pu) — 1 Uk 
with the real part of the pole position. One must sum over k to obtain the number 
density of remnants. 

We shall limit our discussion here to the decay that occurs when the parent ip 
becomes non-relativistic, 

^ fe ~M + ^>T. (61) 
This condition is relevant in interesting cosmological problems of the neutron decay 



at the time of nucleosynthesis and GUT X boson decay at baryogenesis |Tq 



Computation of the temperature dependent part of the occupation number nf 
may proceed by deforming the contour of uo integration, in the same way as in the 
discussion of g{j) above. There are then two types of contribution: One is the 
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pole term that gives the usual Boltzmann suppressed contibution of e /3tJfe . When 
mode-summed, it gives the number density, 

rf3 ^ c -i3(M+k 2 /2M) = ( M2L y>/ 2 e ~ M / T (62) 



(2tt) 3 v 2vr 

This is the familiar Boltzmann suppressed formula. 

The second one is contribution from the continuous complex path that gives the 
power law behavior of temperature dependence. A part of this contribution in the 
region u > y/k 2 + Am? is analytically calculable by using r(u) = + 0[m 2 ] , valid 
for a small daughter mass m. It is 

[A 2 roo /*oo 

n ~ - ; — tttt / dkk 2 dtu — ■= 

64tt 4 M 3 Jo Jk e"/ T - 1 

1 /i 2 T 4 



2880 M 3 



(63) 



This calculation however ignores complicated logarithmic factors in r{uj). 

We numerically computed all terms including the logarithmic factor in r{uj) along 
with 0[m 2 ] corrections. It turns out that the total contribution is ten times larger 
than the analytic result above; in the m — > limit, 

u 2 T A 

n w 3.8 x j^- . (64) 
The main part of this large contribution comes from 

\uj\ < k. With a dimensionless constant introduced by /i — gM, this gives, relative 
to the photon number density, 

We wrote here the numerical value with Gp the weak interaction constant of mass 
dimensions — 2 ( GFm 2 N 10~ 5 ), as if it were relevant to the neutron decay. 

One may estimate the equal time temperature T eq at which the power contribu- 
tion becomes equal to the Boltzmann suppressed number density, to give 

taking as an example \i = 10~ 5 M , the weak interaction strength. This number is 
in an interesting range to affect nucleosynthesis, but we should keep in mind that 
we did not work out the relevant three body decay, n — > p + e + v e . 
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The physical interpretation of the pole term is a conventional one in terms of the 
remnant created by the inverse decay x + X ~ * V 9 ) with all relevant particles on the 
mass shell, hence suppressed kinematically by the Boltzmann factor e _M//T . On the 
other hand, the contribution from the continuous contour can only be interpreted 
as remnant particles far off the mass shell that may exist in thermal equilibrium. 
The usual kinetic approach such as the Boltzmann-like equation is based on the 
rates computed from S-matrix elements on the mass shell and gives the Boltzmann 
suppressed abundance in equilibrium for M ^> T. Our fully quantum mechanical 
approach yields a different result. 

We note that the local friction approximation is equivalent to the pole model 
(with identification of u 2 + 5u> 2 = (Jftzo) 2 , r\ = — 2Qzo ) that ignores the continuum 
integral around the threshold, hence the model fails to describe the off-shell remnant. 

The effect of the off shell remnants seems to play important roles at least in two 



places: in nuclear matter and in the early universe. We shall refer to [19| for some 
recent attempts to estimate medium effects in nucleus. 

We shall mention another application of immediate interest in cosmology; the 
heavy X boson decay for GUT baryogenesis. It has been argued [18] that there 
exists a severe mass bound of order, 

m x > 0[a x m p i] « 10 16 GeV , (67) 

to block the inverse process of the X boson decay so that generation of the baryon 
asymmetry proceeds with sufficient abundance of the parent X particles. The usual 
estimate of the mass bound mentioned above however is based on the on-shell Boltz- 
mann equation. More appropriate formula in this estimate is our remnant number 
density, 

n x ^O[AxlO- 3 ]g 2 x — . (68) 

m x 

(In a more realistic estimate one should consider the X boson decay into quarks 
and leptons. But for an order of magnitude estimate difference in statistics is not 
important.) With the GUT coupling of g 2 x /^ = 1/40, the equal temperature is 
roughly 

M n A 

q w To ' tF^ 1x1 ° ( 69 ) 

eq 

Thus, at temperature of about a tenth of the X mass the Boltzmann suppressed 
formula is replaced by the power formula. The kinematical condition for baryogenesis 
must be reconsidered in view of our off-shell formula. 
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In summary, we extended previous works on the theory of quantum dissipation 
in the linear harmonic environment, with emphasis on exact treatment using the 
reduced density matrix. When applied to the quantum system that decays via 
interaction with the environment, the power law decay is a generic feature towards 
the asymptotic late time limit. The remnant fraction in the asymptotic limit has 
contribution from particles off the mass shell that does not suffer from the Boltzmann 
suppression factor. 
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